Estimating the value of containment strategies in delaying the arrival time of an 
influenza pandemic: A case study of travel restriction and patient isolation 
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With a simple phenomenological metapopulation model, which characterizes the invasion process 
of an influenza pandemic from a source to a subpopulation at risk, we compare the efficiency of 
inter- and intra-population interventions in delaying the arrival of an influenza pandemic. We take 
travel restriction and patient isolation as examples, since in reality they are typical control measures 
implemented at the inter- and intra-population levels, respectively. We find that the intra-population 
interventions, e.g., patient isolation, perform better than the inter-population strategies such as 
travel restriction if the response time is small. However, intra-population strategies are sensitive 
to the increase of the response time, which might be inevitable due to socioeconomic reasons in 
practice and will largely discount the efficiency. 
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I. INTRODUCTION 

During the past decades, extensive efforts have been 
made to investigate the spread of epidemics. Besides 
various epidemiological models having been proposed to 
explore virus transmission in a closed population[lJ, the 
study of network spreading uses structured populations 
to understand the evolution of epidemics in more realis- 
tic social settings These studies have contributed a 
great deal of insightful findings, such as the absence of 
epidemic threshold in scale- free networks 6L the reaction- 
diffusion process, and metapopulation[7, 8|, to name a 
few. These significant advances have raised new issues 
on how to limit or control the spread of infectious dis- 
eases in human society. 

To curb the spatial spread of diseases from city to 
city, a variety of strategies are recommended according to 
World Health Organization(WHO) or United States(US) 
response plans 0: (i) Vaccination of prior groups or dy- 
namic mass vaccination; (ii) antiviral drugs for prophy- 
laxis and treatment; (iii) community-based prevention 
and control; and (iv) travel-related containment mea- 
sures. Except for the fact that travel-related measures 
are implemented at the inter-city level, other strategies 
are mainly performed at the intra-city level. The first two 
pharmaceutical interventions cut down the number of po- 
tential susceptibles or allay the virus transmission rate, 
respectively. Community-based strategies might affect 
individuals (e.g., patient isolation, self-isolation, quaran- 
tine), groups, or entire communities (e.g., cancellation 
of public gatherings, school closures) in a city. Travel- 
related measures mainly result in the restriction or can- 
cellation of nonessential trips. 

By supposing that the outbreak of a pandemic is un- 
derway, many works have studied the efficiency of strate- 
gies by using the metapopulation model, which harnesses 
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the reaction-diffusion framework to sketch human daily 
contacts and mobility. The epidemic reaction takes place 
inside each subpopulation due to personal contacts, and 
the infectious disease cascades subpopulation by subpop- 
ulation via the travel of individuals (here each city is rep- 
resented by a subpopulation). The importance of various 
strategies in decreasing the attack rate or prevalence has 
been extensively studied in Refs. [13, [H| mainly by com- 
putational simulations. Particularly, by analyzing the 
delay of arrival time of the disease |12Hla |. it has been 
shown that the efficiency of travel restriction in slowing 
down the international spread of pandemic influenza is 
limited. 

In these seminal works, the intra- and inter-population 
interventions are seldom compared with each other to 
provide a holistic picture about their value in delaying 
disease invasion. This should give us pause for thought. 
Whether it is reasonable to discard the tactic of travel 
restriction might also depend on how good the intra- 
population strategies perform. In an attempt to study 
this issue, we theoretically analyze the efficiency of two 
kinds of typical containment strategies, namely, travel re- 
striction and patient isolation, which are implemented at 
the inter- and intra-population levels, respectively. We 
mainly use a simple phenomenological model following 
Refs. [ill, Ht| . which considers the importation of an 
infectious disease from a source to a region at risk dur- 
ing the early stage of a pandemic outbreak. Since the 
spreading process cascades subpopulation by subpopula- 
tion, this two-subpopulation version [l5j is a simple model 
but rational approximation of the initial stage of the pan- 
demic. We mainly focus on the impact of strategies to 
delay the arrival time of disease in the subpopulation 
at risk, because no outbreak will occur in an unaffected 
region before the introduction of infectious seeds. Af- 
ter the disease lands in the subpopulation, the ongoing 
endogenous transmission will become the mainstream of 
infections 1 1 1[ 1 1 5l| . Thus the first arrival time of infectious 
travelers is an important quantity characterizing the tim- 
ing of the disease outbreak [ID, [l8l [l9|. 
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II. MODEL DESCRIPTION 

To build the model, we first specify the mechanism of 
individual mobility between subpopulations cc, y. Follow- 
ing Refs. [IM1, at every time step, each individual may 
travel from his current location x(y) to a neighboring sub- 
population y{x) with a per capita diffusion rate uj X y{ui yx )- 
We define the unit time as 1 day. The model proceeds 
with discrete time steps. In reality, the amount of trans- 
portation flows, e.g., air traffic, between cities is often 
symmetric [lj|-[2Cj, which indicates a detailed balance for 
the traffic flows. For simplicity, we assume that the sub- 
populations x, y have the same population size N x — 
N y = N and diffusion rate U) xy = u) yx = uj. Thus there 
are on average uiN individuals departing from each sub- 
population per day. Note that relaxing these two restric- 
tions does not change the main results of this Brief Re- 
port as long as we maintain a detailed balance condition. 
While mobility couples different locations, the epidemic 
reaction process occurs in each subpopulation, where the 
population is mixing homogeneously. We consider a stan- 
dard susceptible-infective-removed (SIR) compartment 
model to represent the influenza-like illness[3, H, [To| . 
At a given time t, the number of susceptible, infec- 
tious, and recovered individuals in x(y) are defined as 
S x (t),I x {t),R x (t)(S v (t),I v (t),Ry(t)), respectively. The 
SIR reaction is governed by the transition rates \x and 
/3[l|. In a unit time, an infectious one recovers and be- 
comes immune at the rate The parameter /3 character- 
izes disease transmissibility, which reflects the combined 
factors of the virus transmission rate and individual con- 
tact rate per unit time[j|. A susceptible individual might 
acquire infection by contact with infectious ones staying 
in the same subpopulation. With the mean-field approxi- 
mation, at time i, the probability for a susceptible one in 
subpopulation x(y) to acquire infection is found by multi- 
plying the density of infectious I x (t)/N(I y (t)/N) by /3[l[. 
In this baseline case, the transfer of susceptible and in- 
fectious individuals is ruled by the diffusion rate uj. The 
epidemic threshold is determined by the basic reproduc- 
tive number Rq = /3/fi, which identifies the expected 
number of secondary infections produced by an infected 
individual during his infectious period in an entire sus- 
ceptible population[l|. 

We next specify the dynamics under interventions. 
Since many socioeconomic factors might defer the im- 
plementation of strategies, we define a response time to 
representing the time interval between the actual incep- 
tion of an outbreak and the time when the strategies 
become available. Travel restriction (TR) mainly affects 
individual mobility between two subpopulations. We de- 
fine the parameter a as the intensity of TR, which means 
that a reduction of fraction a in travel begins at time to, 
i.e., in the model, we decrease the diffusion rate from uj 
to (1 — a)uj after time to- 
Patient isolation (PI) mainly impacts individual com- 
partment transitions. The effect of PI may relate to en- 
forcement by local authorities, or is attributed to the self- 



isolation of infected individuals. For simplicity, we do not 
distinguish between these two aspects. The parameter 77 
is defined to reflect the intensity of PI. It means that on 
average a fraction rj of infectious persons will be isolated 
per unit time after to- We introduce the PI by adding 
an isolation process that each infectious one has a likeli- 
hood to be isolated with rate 77 per unit time. Since these 
isolating individuals have little chance to cause infection, 
we remove them as long as they are isolated. 

III. ANALYTICAL AND SIMULATION 
RESULTS 

Initially, an infectious individual is introduced into 
subpopulation x. Thus the initial condition is 
-^r(0)=l,/y(0)=0. We first analyze the efficiency of TR in 
slowing down disease invasion to subpopulation y. The 
key issue is to evaluate its impact on delaying the first ar- 
rival timc(FAT) of infectious travelers from x. With the 
Poisson process assumption that the diffusion of any indi- 
vidual is independent from that of others, the probability 
that the first infectious individual arrives in subpopula- 
tion y at time t y = t is 

t-i 

p(ty = t) = [1 - (1 - [] (1 <o) ImM , (i) 

ti=l 

which describes that at least one successful transfer of 
infectious individuals from subpopulation xtoy occurs at 
time t, and none at previous time steps (l8l.[l9j. In reality, 
it is general that the number of travelers per day is several 
orders of magnitude smaller than the total population of 
a city, where only small amounts of people leave to travel 
per day. Empirical evidence of worldwide or US domestic 
air transportation^ suggests that the daily diffusion rate 
of individuals on each flight route is of the order 10~ 4 
or less. We here assume uj = 10~ 4 ,iV = 10 6 . Using 
the Taylor expansion, Eq. ([1]) becomes P(t v — t) = 
ul x (t) exp[-u; J2o<u<t T * ■ 

Based on many seminal works [Tol - [l6j |. we assume a 
pandemic influenza with Ro = 1.75 and the infectious 
period = 3 days. In this case, the Malthusian param- 
eter A, the real-time exponential growth rate at the early 
stage of an outbreak[2l],[22|], is /3 — (1 = 0.25. Since w<A, 
the SIR reaction happens at a time scale much faster 
than the diffusion process, thus the number of infectious 
individuals in subpopulation x grows sufficiently before 
subpopulation y is invaded. Meanwhile, at this early 
stage, the infectious ones only make up a small fraction of 
the total population in x,I x (t) <CiV. With a mean-field 
approximation for the evolution of infectious individu- 
als, we haveH H El I X (U) ^ 4(0) exp(A^), U < f. 
Using the continuum approximation Y^,o<t <t ^x{pi) — 
f Q (ItI x (t), we obtain the probability density of FAT, 
P(t) = wexp[Ai — (w/A) exp[AiYL with the mean value 
<t F >= (l/A)(ln(A/w) -7)[ll[l3, where 7 is the Euler 
constant. With the above given parameters, this charac- 
teristic time scale of FAT is <t F >~ 29 days. 
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FIG. 1: (Color online) The analytical results of the relation between the delay of FAT, At, and the intensity of strategies, (a) 
Travel restriction, (b) Patient isolation. The colored squares, circles, and diamonds refer to the cases of to=0, 10, and 20 days, 
respectively. 



In the TR scenario, when the FAT is smaller than 
the response time to, the probability density of FAT 
is still P(t); however, when the FAT is larger than to, 
this probability density becomes P a (t) = (l — (1 — 

(i - «)«)'-(«>) rw* (i - ") IAU) n t0 < tj<t (i - (i - 

a )uy*(h) ~ (1 - a )wexp[Ai - (1 - at)w exp(At)/A - 
au exp(Xto) / X]. We numerically calculate the average 
FAT through < t F >= J Q to tP{t)(1t + TP a (r)dT, and 
get the delay of FAT, At (a), by solving 

At(a) =<t F > - <t F > . (2) 

If the response time t is negligible (to — 0), Eq.fl3]) is 
simplified as 



Ai(a)|to=o = -ln(l-a<)/A, 



(3) 



which recovers the results obtained by the cumulative 
probability P(P < t) in Refs. [lE[l|. Note that Eq © 
is independent from the values of w,N. With A=0.25, 
unless the intensity a is increased to an unpractically 
high level (a > 0.97), At (a) cannot be longer than 2 
weeks. 

To study the PI scenario, we first consider the 
case where the FAT is larger than to- At this early 
stage, we still have the approximation I x (ti) — exp(Ai;) 
when time ti < to; after to, the Malthusian pa- 
rameter becomes X v = A — r\, and thus we have 
I' x (tj) ~ exp(ryto) exp(X v tj) when to < tj < t y . The 
probability density in this case is P v (t) — (1 — (1 — 

^) 4(t) )n < ti <to( 1 - ^ (ti) rw t (i - "Y'^ * 

ul' x (t) cxp(- W / O to I x (r)dr) exp(-c J t [ I' x (t)oIt) 
lo exp[0(to)] exp[A^t — wexp(77to + A^t)/A, ; ], where 
8(to) = r)to — oj exp(Xto)/ X + to cnp(Xto)/ X n . If the re- 
sponse time is negligible (to=0), we simplify the former 
expression as P n (t)\to=o — wexp[A r) t — u exp(X r] t) / X v ], 
which leads to the average FAT, < t F > | to=0 = 
Jo 00 TP v( T )\to=odT ~ (ln(X n /uj) - 7)/A n . In this case, 



we get the relation between At and rj by solving the 
equation 



At(r))\ to = =<t F > | tQ=Q - <t F > 



(4) 



If to > 0, the average FAT is numerically integrated via 
the equation < t F >= J Q to TP(r)dr + /~ rP„(r)dr. We 
therefore have the relation between At and r] as 



At(r)) =<t F > - <t F > 



(5) 



With Eq. Q and A = 0.25, we find that an intermedi- 
ate level of the strategy intensity r\ = 0.12 can adequately 
suspend the arrival of disease to subpopulation y for more 
than 3 weeks. When the response time to = 0, we con- 
clude that the strategy of PI performs better than the 
TR. This is mainly because the TR alone can not miti- 
gate the initial exponential growth of infectious ones in 
the source. However, the strategy of PI is highly sensitive 
to the increase of the response time to- As shown in Fig. 
[TJ when to increases from to 20 days, there is an evi- 
dent decline for the delay At(rj) in the PI scenario, while 
the delay At (a) actualized by implementing the TR is 
robust to the increase of to. 

We further use the dynamic Monte Carlo method to 
simulate the epidemic evolution under different inter- 
ventions. The simulations are performed with discrete 
time steps, and we update each individual's behavior in 
parallel per unit time. The parameters are N = 10 6 , 
uj = 10~ 4 , i?o = 1.75, and yT 1 — 3 days. Initially, an 
infectious individual is introduced into subpopulation x, 
and thus the initial condition is / x (0)=l,/j,(0)=0. When 
the containment strategies are excluded, the epidemic re- 
action and diffusion at each unit time proceed as follows, 
(i) Reaction: Inside each subpopulation, individuals are 
mixing homogeneously. At time t, the probability for any 
susceptible in subpopulation x(y) to acquire infection is 
f3I x (t)/N(f3I y (t)/N). The number of new infections in 
x(y) at time t is extracted from a binomial distribution 
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FIG. 2: (Color online) The relation between the delay of 
FAT, At(a), and the intensity of travel restriction. The gray 
stars are the analytical results with to = 0. The other colored 
symbols are the simulation results with various response times 
t =0, 10, 15,20, and 25 days. 



with probability pi x (t)/N(j3I y (t)/N) and the number of 
trials S x (t)(S y (t)). The number of recovered individuals 
in x(y) is also extracted from a binomial distribution with 
probability /i and the number of trials I x (t)(I y (t)). (ii) 
Diffusion: After all individuals have been updated for the 
reaction, we simulate their diffusion. The number of sus- 
ceptible travelers departing from each subpopulation per 
unit time is also extracted from a binomial distribution 
with probability uj and the number of trials S x {t)(S y (t)). 
The number of infectious and recovered travelers is ob- 
tained in the same way. 

We first study the effects of TR in delaying the arrival 
of disease to subpopulation y. To assemble this factor 
into the simulation, we rescale the per capita diffusion 
rate w by a multiplier 1 — a, where the parameter a re- 
flects the intensity of TR. The strategy is activated after 
a given response time to- FigureOprovides a holistic view 
about the relation between the delay of FAT, At(a), and 
the restriction intensity a. Since the disease might die 
out due to randomness, every data point is obtained by 
averaging the simulations with the successful transfer of 
infectious ones among 10 4 times of Monte Carlo random 
experiments, each of which is simulated with 500 time 
steps. The gray stars are the analytical results obtained 
by Eq. @, which agree well with the simulations. When 
a = 0.3,0.6, and 0.9 and to = 0, the simulations show 
that At(a) ~ 2,4, and 11 days, respectively. Even if the 
restriction intensity is elevated to an unpractically high 
level, e.g., a — 0.97, At(a) is still less than 3 weeks. It 
is clear that At is small if the time scale of the initial 
exponential growth 1/A is small [see Eq. [3]. We further 
study the impact of the response time on the efficiency 
of TR. In Fig. [21 unless to approaches < t F >, which is 
the average FAT without TR, and a is large, there is no 
evident decline for the simulation results of At (a). 

We next study the effects of PI in delaying disease in- 
vasion. To introduce this factor in the model, we add 



an isolation process before the reaction process at each 
time step after to. The parameter r\ reflects the inten- 
sity of PI. Per unit time, the number of newly isolated 
individuals in subpopulation x(y) is extracted from a bi- 
nomial distribution with probability 77 and the number of 
trials I x (t)(I y (t)). Figure |3J a) presents the relation be- 
tween the delay of FAT, At, and the isolation intensity 77 
with to = 0. For each 77, we perform 10 4 times of Monte 
Carlo random experiments, each of which is simulated 
with 500 time steps. Due to the randomness embedded 
in the dynamical process, the infectious individuals in 
source x might be totally eradicated before traveling to 
subpopulation y. With a given 77, we measure Atiif) by 
averaging the simulations that the infectious ones from 
source x successfully jump to subpopulation y. The re- 
sults are highlighted by the red squares in Fig. G2a). The 
gray stars are the analytical results obtained by Eq.Q. 
If the isolation intensity 77 is at a small or intermediate 
level(77 < 0.18), the simulation results agree well with 
the theoretical predications. However, if the intensity 
77 is extremely large, the simulations obviously deviate 
from the analytical results. In this latter case, since the 
Malthusian parameter A,, is quite small, there is a huge 
likelihood of totally eradicating the infectious individuals 
at the early stage of an outbreak due to randomness. For 
instance, when 77 = 0.2, 0.22, the fraction of eradication in 
all independent modeling realizations reaches 97.7% and 
99.2%, respectively, while for 77 = 0.12, the fraction of 
eradication is only 75.6% [see the dark cyan diamonds in 
Figj3ja)]. If 77 > 0.25, the Malthusian parameter A < 0, 
the disease hardly persists in the population. With the 
same condition that to = 0, the strategy of PI is more 
efficient than TR: An intermediate level of isolation in- 
tensity 77 can adequately delay the arrival of disease for 
about 1 month. 

Figure [3jt>) shows the impact of the response time to 
on the delaying effects of PI. For a small t , e.g., t = 10, 
which is much smaller than <t F >, an intermediate level 
of PI(e.g., 77 = 0.14) still suspends the arrival of disease 
for about 3 weeks. This achievement exceeds the per- 
formance of TR even with an extremely high restriction 
intensity. The simulations also illuminate that the PI 
is sensitive to the increase of to- There is a remarkable 
decline in the simulation results of At (77) when to ap- 
proaches < t F >. For instance, 25 days of waiting to 
implement the strategy (to = 25) will only postpone the 
arrival of disease in subpopulation y for about 2 weeks 
at most. 

Actually, other intra-population interventions can also 
be analyzed under this framework. For instance, social 
distancing limits public activities to reduce personal con- 
tacts, which can be reflected by rescaling the disease 
transmission rate /3 with a multiplier 1 — tp when time 
t > to- At the initial stage of an outbreak, the Malthu- 
sian parameter becomes X v = (l — tp)/3 — fi. From a math- 
ematical point of view, we can adjust the parameters tp, 77 
to allow A^ = A^. Therefore, the above analysis can cover 
this scenario. 
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FIG. 3: [Color online) The effects of patient isolation in delaying disease invasion, (a) The relation between the delay of 
FAT, At (77), and the isolation intensity n with to = 0. The dark cyan diamonds show the fraction of eradication. (b) The 
simulation results with to = 0, 10, 15, 20, and 25. 



IV. SUMMARY 

In sum, the intra-population interventions, e.g., pa- 
tient isolation, perform better than the inter-population 
strategies such as travel restriction if the response time 
is small. Therefore, the intra-population strategies are 
more beneficial in delaying the spatial spread of pan- 
demic influenza if they are implemented very promptly. 
However, the intra-population measures are sensitive to 
the increase of response time, which might be inevitable 
due to miscellaneous socioeconomic reasons in reality and 



largely discounts the efficiency. 
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